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Abstract 

We present the design, implementation and data taking performance of the Microwave Detection of Air 
Showers (MIDAS) experiment, a large field of view imaging telescope designed to detect microwave radiation 
from extensive air showers induced by ultra-high energy cosmic rays. This novel technique may bring a 
tenfold increase in detector duty cycle when compared to the standard fluorescence technique based on 
detection of ultraviolet photons. The MIDAS telescope consists of a 4.5 m diameter dish with a 53-pixel 
receiver camera, instrumented with feed horns operating in the commercial extended C-Band (3.4 - 4.2 GHz). 
A self-trigger capability is implemented in the digital electronics. The main objectives of this first prototype 
of the MIDAS telescope - to validate the telescope design, and to demonstrate a large detector duty cycle - 
were successfully accomplished in a dedicated data taking run at the University of Chicago campus prior to 
installation at the Pierre Auger Observatory. 



1. Introduction is 

19 

The origin and composition of Ultra-High En- 20 

ergy Cosmic Rays (UHECRs) remains uncertain pQ , 21 

even after the progress made by the latest genera- 22 

tion of experiments [2J [3] . Due to the strong flux 23 

suppression above 10 19 eV [H [5], very large de- 24 

tection areas are necessary to study cosmic rays 25 

at these energies. A future UHECR Observatory 2 e 

based on standard techniques - Surface Detector ell"- 27 

rays (SD) and Fluorescence Detectors (FD) - may 28 

be limited by cost and difficulty of deployment. In 29 

this context, radio detection techniques are attrac- 30 

tive thanks to the low cost of individual elements, 31 

the little maintenance required and a nearly 100% 32 

detection duty cycle. 33 

Radio emission in the MHz range from extensive 34 

air showers (EAS) has been actively studied in the 35 



last decade by the LOPES [5] and CODALEMA [7] 
experiments. MHz radio-detection is now well es- 
tablished, and AERA [5], in commissioning phase 
at the Pierre Auger Observatory, will have a suffi- 
ciently large instrumented area to explore the po- 
tential of this technique for detection of the highest 
energy cosmic rays. 

The use of the microwave (GHz) band for EAS 
detection was originally pursued by Jelley, Char- 
man and collaborators [3 [TU] in the late 1960s, but 
abandoned due to the lack of satisfactory under- 
standing of the emission mechanisms and to lim- 
itations in detector technology. Thereafter it re- 
mained mostly unexplored, until recent laboratory 
measurements [11] with particle beams have re- 
newed the interest in this part of the radio spec- 
trum. These measurements suggest that microwave 
radiation is emitted from the weakly ionized air 
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Figure 1: The MIDAS telescope at the University of Chicago, 
with the 53-pixel camera at the prime focus of the 4.5 m 
diameter parabolic dish reflector. 

36 plasma of free electrons produced by the EAS in- 

37 duced ionization of the atmosphere. The radia- 

38 tion is expected to be continuous and relatively flat 

39 in frequency unpolarized and emitted isotropically, 

40 and its intensity to scale with the number of parti- 

41 cles of the shower. 

42 Detection of an isotropic emission in the GHz 

43 range - akin to the detection of ultraviolet fluo- 

44 rescence photons by FD pioneered by the Fly's 

45 Eye experiment |12j and currently used by the 

46 Pierre Auger Observatory [T3] and the Telescope 

47 Array [T3] - allows for the measurement of the EAS 

48 development in the atmosphere, which provides a 

49 calorimetric measurement of the energy and crucial 
so information on the mass composition of primary 

51 cosmic rays. A GHz radio telescope would over- 

52 come the limitations of the FD technique, i.e. data 

53 taking only during moonless nights (sa 15% duty 

54 cycle) and significant systematic uncertainties in- 

55 troduced by light attenuation in the atmosphere. In 

56 fact, microwave detectors can operate 100% of the 

57 time, and attenuation in the GHz range is minimal, 

58 even with rain or clouds. Moreover, commercial 

59 off-the-shelf GHz equipment, mostly developed for 

60 satellite TV reception, is readily available and in- 

61 expensive. Microwave telescopes could provide ex- 




Figure 2: A front view of the receiver camera, with LNBFs 
closely packed to maximize the sensitivity over the focal 
plane. 

62 isting UHECR experiments of unprecedented sen- 

63 sitivity to primary comic ray composition, and be 

64 employed in a future large scale observatory. 

65 Several complementary approaches to microwave 

66 detection of EAS are currently being pursued, in- 

67 eluding the AMBER and EASIER detectors Q2] at 

68 the Pierre Auger Observatory, and the CROME ex- 

69 periment at KASCADE [16 . Also, new laboratory 

70 measurements with particle beams are being per- 

71 formed [T71IIH] to better characterize the microwave 

72 emission. 

73 In this paper, we present the Microwave Detec- 

74 tion of Air Showers (MIDAS) experiment, an imag- 

75 ing telescope whose primary objective is to confirm 

76 the microwave emission from EAS, and to demon- 

77 strate the feasibility of a low cost design for this 

78 novel technique. The MIDAS telescope is described 

79 in Sec. 2. The electromagnetic simulations of the 
so telescope response are illustrated in Sec. 3. The 
si calibration procedures and the measured sensitiv- 

82 ity of the instrument are presented in Sec. 4. A 

83 realistic simulation of the MIDAS detection of EAS 

84 is described in Sec. 5. Operation and data taking 
as performance of the MIDAS telescope are presented 
ae in Sec. 6, and conclusions are drawn in Sec. 7. 

87 2. The MIDAS telescope 

as The telescope consists of a large parabolic dish 

89 reflector with a receiver camera at its prime focus, 

90 installed on the roof of the Kersten Physics Teach- 

91 ing Center at the University of Chicago (Fig. [T]) . In 

92 the following, the design and technical implemen- 

93 tation of the different components of the MIDAS 

94 telescope are described. 
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Figure 3: The analog electronics chain. See text for details 
on the different components. 



2.1. Reflector and receiver camera 

The parabolic dish reflector (Andrew) has 4.5 m 
diameter and f/D — 0.34. A motorized alt-azimuth 
mount allows for telescope movements in a range of 
90° in elevation and 120° in azimuth. The remote 
control of the telescope pointing (0.1° precision) is 
integrated in the data acquisition (DAQ) system 
(Sec. 2.5 ). A 53-pixel receiver camera is mounted on 
the prime focus of the dish, covering a field of view 
of about 20° x 10°. The microwave receivers are 
arranged in seven rows, and staggered to maximize 
the sensitivity across the focal plane (Fig. [2J . 

A commercial low noise block feed horn (LNBF) 
operating in the extended C-band (3.4 - 4.2 GHz) 
is used for the receiver. These feeds (WS Inter- 
national) are mass-produced for consumer satellite 
television. The LNBF integrates a feed horn, low 
noise amplifiers, and a frequency downconverter. 
The feed can receive two orthogonal linear polar- 
izations which are remotely selectable through the 
LNBF power voltage level setting. A 5150 MHz lo- 
cal oscillator in the frequency downconverter mixes 
the input RF signal down to a frequency interval of 
950-1750 MHz, which is transmitted with minimal 
loss through standard coaxial cable. The receiver 
bandwidth, its gain T and noise temperature were 
measured to be about 1 GHz, 65 dB and 20 K, re- 
spectively. 

Power to the LNBF is provided through 30 me- 
ters of commercial quad-shielded RG-6 coaxial ca- 
ble, which also brings the RF signal, after ampli- 
fication and downconversion, from the telescope to 
the counting room. 

2.2. Analog electronics 

The analog electronics chain is summarized in 
Fig. [3j The RF signal in the coaxial cable is first 



Figure 4: An analog electronics tray with eight channels. 
The bias tee (f), impedance adapter (2), power detector (3), 
and DC voltage inputs (4) are visible. 



passed through a bandpass filter (1.05-1.75 GHz). 
The purpose of the filter is to reject interference 
from radar altimeters of airplanes, which was iden- 
tified as a major source of background in the early 
stage of commissioning of the telescope. A power 
inserter (commonly called bias tee) provides the 
DC voltage to the LNBF, while letting through 
the RF signal. An impedance adapter matches the 
standard 75 Ohm impedance of these commercial 
satellite TV components to the 50 Ohm impedance 
of the RF power detector (Mini-Circuits ZX47-60- 
S+). The overall power loss L in the 30 m of cable 
and in the analog electronics was measured to be 
17 dB. 

The power detector (input bandwidth w 8 GHz) 
responds logarithmically to an RF power P in the 
range -55 dBm to 5 dBm, with a voltage output 
between 2.0 V and 0.5 V: 



V = Vb-lOalog(P), 



(1) 



with P in mW. The characteristics of the 53 power 
detectors were individually measured, with typical 
values around V a = 0.625 V and a = 0.025 V/dB. 
Also, their response time was found to be about 
100 ns, well suited for typical pulses of microseconds 
duration expected from an EAS crossing the field 
of view of a pixel. 

The analog electronics components are arranged 
in trays - 8 channels per tray - for distribution of 
DC power and routing of the signals to the digital 
electronics (Fig. |4| . 

2.3. Front-End electronics and digitizers 

The signal from the power detector is digi- 
tized with 14-bit resolution (calibration constant 




Figure 5: Digital electronics board with 16 Flash-ADC chan- 
nels. The module includes an FPGA for the First Level Trig- 
ger logic and a VME interface. The board was designed by 
the Electronics Design Group at the Enrico Fermi Institute 
of the University of Chicago. 



164 b = 7 ADC/mV) at a sampling rate of 20 MHz by 

165 custom-made Flash-ADC (FADC) boards (Fig. p)}. 

166 Up to 2048 samples are stored in a circular buffer 

167 and processed by a first level trigger algorithm 
lea implemented in the on-board FPGA. A standard 
leg VME interface allows for board control and data 

170 readout. 

171 Four FADC boards, each with 16 channels, are 

172 sufficient for the entire receiver camera. A Master 

173 Trigger Board, also equipped with FPGA logic and 

174 VME interface, provides the global clock for the 

175 FADC synchronization and performs a high level 

176 trigger decision. Each FADC board is connected to 

177 the Master Trigger board through Low Voltage Dif- 

178 fcrential Signaling (LVDS) lines carrying the clock 

179 and trigger signals. 

180 A VME module equipped with a GPS receiver 

181 (Hytec 2092) tags the time of the triggered events 

182 with 10 ns precision. Details of the MIDAS digital 

183 electronics can be found in 1191 . 
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Figure 6: Illustration of the FLT. The digitized time 
trace for a 5 fis RF pulse from the calibration antenna (see 
Sec. |2,5| , with the ADC running average of 20 consecutive 
time samples superimposed as a gray histogram (red in the 
color version). An FLT is issued when the running average 
falls below the threshold, indicated by the horizontal line. 
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Figure 7: Illustration of the SLT. The 4 highlighted pixels 
(top panel), located in three different FADC boards, have 
an FLT occurring in the same 50 ns time sample. Each 
FADC board sends through the LVDS lines the FLT status 
of its pixels to the Master Trigger Board (bottom panel, left), 
where a matching SLT pattern is found. An SLT is then sent 
back through the LVDS lines to the FADC boards (bottom 
panel, right) triggering the entire camera for data readout. 
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Figure 8: Monitoring data from one pixel for a six-hour data taking run. On the left, the ADC baseline averaged over 10 ms 
(upper panel) and the FLT threshold (bottom panel) as a function of time. The threshold follows closely the variations of the 
baseline, keeping the FLT rate during the run (right panel) around 100 Hz. 



184 2.4- Trigger 21s 

219 

185 MIDAS implements a multi-level trigger architec- 

186 ture, optimized for transient events with topology 

187 and time structure compatible with signals from an ^ 

188 EAS. The system is modeled on the one success- 

J 223 

189 fully operating at the FD of the Pierre Auger Ob- 

190 servatory [13] . For each channel, a First Level Trig- 225 

191 ger (FLT) algorithm identifies pulses in the FADC 

192 trace. The FLT information from all channels is 

193 sent to the Master Trigger Board where the Second 226 

194 Level Trigger (SLT) searches for patterns of FLTs 227 

195 compatible with a cosmic ray track. 22a 

196 The FLT status bit of a given channel is acti- 229 

197 vated whenever the ADC running sum of 20 consec- 230 

198 utive time samples falls below a threshold (Fig. [6]). 231 

199 To compensate for changes in the noise level, the 232 

200 threshold is continuously regulated in order to keep 233 

201 a stable FLT rate around 100 Hz. The active status 234 

202 of the FLT bit is extended for 10 /us after the trig- 235 

203 ger, to allow for time coincidences between different 236 

204 channels in the Master Trigger Board. 237 

205 Every 50 ns, each FADC board transmits 16 bits, 23a 

206 corresponding to the FLT status of its 16 channels, 239 

207 to the Master Trigger Board (Fig. [7]), where the 240 

208 SLT trigger algorithm searches for 4-fold patterns 241 

209 of channels with FLT triggers overlapping in time. 242 

210 There are 767 patterns compatible with the topol- 243 

211 ogy of a cosmic ray shower (straight tracks across 244 

212 the camera). When an SLT matching pattern is 245 

213 found, an event trigger signal is distributed back to 246 

214 the FADC boards and to the GPS module. At this 247 

215 point, a block of 100 fis of ADC data (including 500 248 

216 pre-trigger samples) is frozen in the memory buffer 249 

217 and made available for readout via VME. All 53 250 



channels, even those not participating to the trig- 
ger decision, are readout. 

Higher levels of trigger control have also been im- 
plemented. For example, the FLT threshold regu- 
lation and data acquisition are inhibited whenever 
the SLT rate is higher than a certain limit, to be au- 
tomatically restored when the SLT rate comes back 
to a normal level. 

2.5. Data acquisition, monitoring and operation 

A single board computer (GE Intelligent Plat- 
forms V7865) acting as VME master is used for 
the data acquisition. The DAQ software monitors 
the relevant VME registers of the FADC boards, 
and reads out through the VME bus the GPS time 
stamp and the data available in the buffer after an 
event trigger. The event data are then assembled 
and written to disk in ROOT format 20 . 

For monitoring purposes, the ADC baseline av- 
eraged over 10 ms (calculated in the FPGA of the 
FADC board), the FLT rate, and the FLT thresh- 
old of each channel are readout and recorded every 
second. An example of monitoring data collected 
during a six-hour data taking run is presented in 
Fig. [8] 

Additional monitoring information is obtained 
from RF pulses periodically illuminating the re- 
ceiver camera. For this purpose, a patch antenna 
with a wide beam (HPBW=70°) is mounted at the 
center of the reflector dish. The antenna is driven 
by an RF signal generator located in the counting 
room, which every 15 minutes produces a train of 
ten RF pulses of 10 fis width. These data provide 
monitoring information for all of the channels at 
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251 the same time, ensuring they are properly operat- 297 

252 ing and that the trigger is working as designed. 29a 

253 Data runs last six hours and are automatically 299 

254 restarted, with data files backed up on a server for 300 

255 offline analysis. A fully automatic and remote oper- 301 

256 ation of the telescope was achieved through a slow 302 

257 control software for the antenna positioning, and 303 

258 for the status of power supplies and of the VME 304 

259 crate. 305 



260 3. EM simulation of telescope efficiency 30s 

309 

261 The efficiency of the MIDAS telescope has a sig- 310 

262 nificant angular dependence, due to the large field 311 

263 of view implemented in a parabolic reflector. For a 312 

264 proper calibration of the telescope and for a realistic 313 

265 estimate of its sensitivity to EAS, the power pat- 314 

266 tern of each pixel - i.e the pixel detection efficiency 315 

267 as a function of the direction (9, (f>) of the incident 316 

268 microwave radiation with respect to the telescope 317 

269 boresight - must be determined. A full electromag- 318 

270 netic (EM) simulation of the reflector and receiver 319 

271 camera has been developed for this purpose. 320 

272 At radio and microwave wavelengths, it is conve- 321 

273 nient to simulate the antenna system (reflector and 322 

274 feed) in emission mode, making use of the so-called 323 

275 Reciprocity Theorem [2TJ to obtain its effective 324 

276 area Aj in a given direction: 325 



At 



\ 2 

= 7- G * 



(2) 



278 where A is the wavelength of the radiation and 330 

279 Gi(9 1 cf>) is the gain of the antenna system for pixel 331 

280 i, defined as the power emitted per unit solid an- 332 

281 gle in the direction (#, <f>) normalized to the corre- 333 

282 sponding power emitted by an isotropic radiator. 334 

283 In practical terms, the antenna effective area sum- 

284 marizes the dependence of the telescope gain on the 

285 dimensions of the reflector, the frequency of the ra- 

286 diation, and overall illumination (i.e. taper) of the 

287 feed. 

288 The feed is simulated as two cylindrical waveg- 

289 uides - the one corresponding to the feed aperture 33a 

290 having a diameter d = 6.7 cm and the other hav- 339 

291 ing a diameter d = 5.5 cm - joined by a small coni- 340 

292 cal throat. The dominant propagating mode is the 341 

293 TE11 (electric field perpendicular to the feed axis) 342 

294 with a cutoff wavelength of 9.4 cm determined by 343 

295 the diameter of the end section (\j*e = 1.71c?) 344 

296 [2 3) . The calculated radiation pattern of the feed in 345 



polar coordinates for the E-plane (</> = tt/2) and H- 
plane ((f) — 0), and the direction of the electric vec- 
tor field over the feed aperture are shown in Fig. [9j 

To obtain the antenna effective area, the electro- 
magnetic field distribution over the telescope aper- 
ture is calculated. For each feed, the telescope aper- 
ture is found by intersection of rays emitted by the 
feed with the plane perpendicular to the telescope 
axis and containing the dish focal point. The paths 
of rays from the feed to the aperture plane are de- 
termined in the optical limit using Fermat's princi- 
ple of least transit time. Spherical wave propaga- 
tion is assumed from feed to dish and plane wave 
propagation is assumed from dish to aperture. A 
proper treatment of the optical paths is essential for 
MIDAS, since the lateral displacement of the feeds 
with respect to the dish focal point produces sig- 
nificant differences in their gain. Also, the shadow 
of the receiver camera, which increases the size of 
sidelobes, is included in the simulation. The radia- 
tion field is then calculated in the far-field limit by 
the Fourier transform of the aperture field distribu- 
tion, and the power is obtained as the square of the 
associated Poynting vector. 

Equation [2] holds for an antenna system com- 
pletely matched to its transmission line. In prac- 
tice, reflections at the feed entrance, at the throat 
between the cylindrical sections and at the coupling 
between the feed and the coaxial line reduce the 
antenna effective area. Their effect can be calcu- 
lated introducing an equivalent transmission line 
model with impedance mismatches which produce 
reflected power waves. Additional signal losses in 
the coaxial cable and in the analog electronics are 
included in the simulation. Also, the reflector does 
not collect all the power emitted by the feed, and 
part of it spills over, further reducing the antenna 
total effective area. 

The result of the EM simulation can be expressed 
as a relative power pattern for each pixel i: 



g<(M) 

Gi 5 (0,0) : 



(3) 



where Gi5(0,0) is the calculated gain of the cen- 
tral pixel in the direction of boresight. The pat- 
tern for the central feed (Fig. 10 left) is symmetric 
around the boresight with relatively small sidelobes. 



Feeds away from the focal point (Fig. 10 right) have 
smaller gains and bigger, asymmetrical lobes. From 
Eq. [2j the effect ive area of the antenna system to 
a microwave flux incident from a given direction is 



G 




Figure 9: EM calculation of the feed for the dominant TEn mode at 3.8 GHz. Left: Normalized electric field radiation pattern 
for the MIDAS feed in polar coordinates. Both the E-plane (<fi = tt/2) and H-plane are shown (cj> = 0). Right: Electric field 
vector over the feed aperture. 




Figure 10: Left: Power radiation pattern for the central feed (pixel 15) of the MIDAS camera. This feed presents the maximum 
antenna gain of the receiver camera, and a symmetrical pattern with very small sidelobes. Right: Radiation pattern for a feed 
laterally displaced by 16.5 cm (~ 2A) from the camera focus (pixel 13). Its gain is ~50% smaller, and presents a significant 
coma lobe. 



346 then given by: 

A i (9,4,)=A eff e i (0,<t>), (4) 

348 where A e ff = j^Gi5(0,0) is the effective area of 

349 the dish. 

350 The overall sensitivity of the camera can be 

351 calculated summing over the pixels: e{9, <fi) = 

352 Yli e i{@ i > an d is presented in Fig. [TT| Even with 

353 the compact arrangement of the MIDAS feeds over 

354 the focal plane, there is a region of reduced effi- 

355 ciency between adjacent feeds. Also, the efficiency 

356 at the borders of the field of view is as low as 20% of 

357 the efficiency at the center. While not uniform over 



358 the camera, the telescope efficiency is reasonably 

359 high in a large portion of its field of view, which is 

360 sufficient for the goals of this first low-cost proto- 

361 type for the MIDAS concept. 

362 4. Telescope Calibration 

363 Dedicated calibration measurements were per- 

364 formed during the commissioning of the telescope, 

365 which were followed up during data taking to verify 

366 the stability of the system. 

367 The Sun, whose flux in the microwave range is 

368 measured daily by several solar radio observato- 
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Figure 11: Radiation pattern of the MIDAS camera, calcu- 
lated by summing over all pixels, as function of the angle 
with respect to the antenna boresight. 



ries around the world, was used as the primary 
source for the absolute calibration of the MIDAS 
telescope. This calibration was cross-checked with 
other known astrophysical sources. 



373 The monitoring patch antenna (Sec. 2.5 ) was also 

374 used for relative channel calibration. 

375 4-1- Absolute calibration 

376 An unpolarized microwave flux density, F, inci- 

377 dent on the MIDAS reflector from a direction (0, <ft) 

378 produces, at the input of the power detector of a 

379 given pixel i, a power Pj given by: 



410 
411 
412 
413 
414 
415 
416 
417 
418 



Pi 



eff 



■ F 



■F, 



(5) 



where ei(9,(j)) is the pixel relative power pattern, 
A e ff is the effective area of the dish, Bi is the 
bandwidth of the electronics (essentially given by 
the pixel bandpass filter), I\ is the gain of the re- 
ceiver, Li is the loss due to the 30 m cable and the 
analog electronics, and the factor 2 takes into ac- 
count the fact that the linearly polarized receivers 
detect only half of the available power. The output 
voltage of the power detector (Eq. [T]) is digitally 
sampled in the FADC board, yielding a signal rii in 
ADC counts: 



= nj-10fcilog [otj • 6i(6, (/)) ■ F] 



(6) 



where n" = bV , and h = ab ~ 175 ADC/dB with 
a, Vo and b as defined in Sec. |2.2| and |2.3[ respec- 
tively. The calibration constant fcj is independently 



419 
420 
421 
422 
423 
424 
425 
426 
427 
428 
429 
430 
431 
432 
433 
434 
435 
436 
437 
438 



measured for each channel (see Sec. 4.2) 



Even when the telescope is pointed towards a re- 
gion of the sky with no microwave sources, a power 



P^ ys is present at the input of the power detector, 
with a corresponding signal in ADC counts. 
Sources for this signal include radiation from the 
sky in the main beam of the antenna, radiation 
from ground in the antenna sidelobes, and electron- 
ics noise, dominated by the receiver noise temper- 
ature. An equivalent flux density, F* ys , is defined 
from P? vs = aiF^ vs . 

A measured signal rii can be converted into an ab- 
solutely calibrated flux density F by deriving from 
Eq.© 



sys 



fOfc.log I 



• F 



Ff 



(7) 



where ej(#,0) is taken from the simulation de- 
scribed in Sec. [3j and n s i ys and the corresponding 
F? must be determined through calibration pro- 
cedures. 

The Sun was used as the main calibration source 
of the MIDAS telescope, by measuring the pixel 
signal during a transit of the Sun in its field of view. 
From Eq. [7J the time evolution of the signal is given 
by: 



rii(t) = n sy: 

10 h log ( I 



E % {6(t)^{t))-F s 



(8) 



where 



.s IIS 
1 M 



is the ADC baseline measured before the 



Sun transit, (6(t), 4>(t)) is the sun position in the sky 
at time t, and the Sun flux density, Fs U n, is taken 
from the daily measurements by the Nobeyama ob- 
servatory 24 . The value of F^ ys which best de- 
scribe the calibration data is then obtained. An 
extensive measurement campaign was performed, 
with each Sun transit aiming to the calibration of 
a row of pixels. An example of a Sun calibra- 
tion run for the top row of pixels is presented in 
Fig. |I2| These measurements not only provide an 
absolute calibration of each pixel through F^ ys , but 
also demonstrate the quality of the simulation of 
ti(0, </>). In fact, simulations were found to describe 
very well the pixel power pattern in both the width 
and the relative efficiency of the lobes, even for pix- 
els at the edges of the camera where strong aberra- 
tions are present. 

As an example, several calibration measurements 
of the central pixel yielded a F° ys = 1.96 • I0 4 J3Q 

*1 Jy = Id" 26 W/m 2 /Hz 
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Figure 12: Example of a calibration measurement of the MIDAS telescope with a Sun transit. Top right inset: Sun trajectory 
in the top row of the camera pixels (azimuth measured clockwise from North). The measured ADC signal (closed dots) for 
each pixel is shown in the other panels, with the superimposed line representing the EM simulation for the best value of F^ ys 
(sec text). 



equivalent, in the common notation of radio astron- 
omy, to a system temperature T sys of: 

Aeff 



T, 



*15 



sys 



2kh 



65 K, 



(9) 



where fc& is the Boltzmann constant, and A e ff — 
9.1 m 2 is the effective area (« 60% of the geometric 462 
area) of the telescope as determined by simulations 
described in Sec. [3] The other pixels had a similar 463 
system temperature. The estimated systematic un- 
certainty on the measurement of F? ys is 15%, dom- 465 
inated by the uncertainty of the Nobeyama mea- 466 
surement and of the telescope pointing. 

Measurements of the Moon (Fmoou — Fsun/WO) " m 
and of the Crab Nebula (F Cra b — -^Sun/1000) were 469 
also performed, providing a cross-check of the ab- 
solute calibration. Using the Sun calibration, the 470 
Moon temperature was found to be 245 K, to be 
compared with 268 K |25j (the measurement was 471 
taken a day after the full moon). The flux density 472 



of the Crab Nebula was measured to be 760 Jy, to be 
compared with 718±43 Jy at 3.38 GHz [26]. Both of 
these results agree to better than 10%, within the 
estimated systematic uncertainty of the Sun cali- 
bration. 

4-2. Pixel calibration constants and timing 

The calibration constant ki in Eq. [6] was mea- 
sured by using the patch antenna mounted at the 
center of the dish. The pulse power P pu ise of the 
RF signal generator driving the patch antenna was 
changed over a wide range, and the corresponding 
signal n p i ulse was measured for each channel. From 
Eq. [7J one expects: 



pulse sys 



10 ki log 1 + 



pulse 



p 



(10) 



where is the signal measured when the RF sig- 
nal generator is off, P? ys is the corresponding power 



9 



12000 - 




173 ± 1 ADC/dB 



7000 



482 
483 
484 
485 



-40 



-30 



-20 -10 

P p u, S e(dBm) 



10 



20 



490 
491 
492 
493 
494 
495 



Figure 13: Example of a measurement of the calibration 
constant k. Closed dots represent the measured signal in 
pixel n. 32 for different powers of the calibration antenna 
pulses. The line is the result of a fit to derive k%2- For 
Ppulse greater than 8 dBm, measurements are affected by 
saturation of the LNBFs and were excluded from the fit. 
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Figure 14: Distribution of the measured calibration con- 
stants k. 



at the input of the power detector, and fi ■ P pu ise is 521 

the power induced by the calibration pulse at the 522 

input of the power detector. The factor fi may 523 

change from pixel to pixel, depending on the dis- 524 

tance of the pixel from the patch antenna, on the 525 

relative orientation of the linear polarizations of the 526 

emitter antenna and the receiver, and on the signal 527 

loss L. An example of a fit of Eq. [lO] , with fej, 52a 

n^ vs and fi/Pi VS as free parameters, to the calibra- 529 



tion data for one channel is shown in Fig. 13 The 
spread of the distribution of the fitted calibration 



constants for the 53 channels (Fig. 14) is 2.6%, with 
an average of k = 172 ADC/dB. 

The calibration pulses have a fast risetime 
(<1 ns) and illuminate all the camera receivers si- 
multaneously, allowing for a measurement of their 
time response and synchronization. The time at 
which the amplitude of the measured pulse reaches 
50% of its maximum value, t$o, was taken as an esti- 
mator. The distribution of £50 was found to have an 
RMS value of 25 ns, more than adequate for the typ- 
ical pulses of microseconds duration expected from 
an EAS crossing the field of view of a pixel. 

5. Sensitivity to EAS 

The sensitivity of the MIDAS telescope to mi- 
crowave emission from EAS has been studied with 
Monte Carlo simulations. The shower development 
in the atmosphere - i.e. the number of charged par- 
ticles N at any given atmospheric depth - is simu- 
lated with a Gaisser-Hillas parameterization [27] . 
including proper fluctuations of the shower pro- 
file. The shower arrival direction is isotropically 
distributed, with its landing point randomly dis- 
tributed over a large area around the telescope. For 
a given shower geometry, the number of particles, 
N(t), along the shower profile is calculated, where 
t is the time in 50 ns samples of the MIDAS digital 
electronics. 

The microwave flux density at the detector aper- 
ture is then calculated, following the model in [IT] , 
as: 



F(t) = F Icf 



Po 



mi ) J ' (jvj) ' 



where F re { is the flux density at a distance d = 0.5 m 
from a reference shower of E re f — 3.36 x 10 17 eV, 
R(t) is the distance between the detector and the 
EAS segment, and p(t) (po) is the atmospheric den- 
sity at the altitude of the EAS segment (at sea 
level). -/V rc f is the average number of shower parti- 
cles at the maximum of the EAS development for 
a proton primary of energy E re f. The parame- 
ter a characterizes phenomenologically the coher- 
ence scaling relationship for the EAS microwave 
emission, with a = 2 (a = 1) corresponding to a 
fully coherent (incoherent) emission. Laboratory 
measurements [11] suggest a reference flux density 
F f ° rcf = 1.85 x 10" 15 W/m 2 /Hz and full coherence 
at shower maximum. 
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Figure 15: Event display of a 3 • 10 19 eV simulated shower 
landing approximately 10 km from the telescope. In the top 
panel, pixels with an FLT are highlighted, with color coded 
by arrival time. In the bottom panel, ADC running aver- 
ages of 20 consecutive time samples for the selected pixels 
(identified by black dots in the top panel) arc shown. The 
running average of each pixel is referred to the threshold 
level (horizontal line) for display purposes. 

530 In order to convert the microwave flux density at 

531 the detector aperture into a signal in ADC counts, 

532 the efficiency maps and calibration constants de- 

533 scribed in Sec. [3] and Sec. [4] were implemented in 

534 the simulation. For each channel, the actual value 

535 of n sys and its fluctuation were taken to be equal 

536 to their average values measured during several 

537 months of data taking, providing a realistic sim- 

538 ulation of the telescope sensitivity. The FLT and 

539 SLT algorithms of Sec. |2.4| were also implemented, 

540 and all simulated events fulfilling the SLT condition 

541 are written to disk in the same format as the data. 

542 An example of event simulated with F rc [ = F® ei and 

543 a — 1 is shown in Fig. |15| 

544 Simulations with different assumptions on the 

545 characteristics of the microwave emission from EAS 

546 were performed. For F TC { — F® e[ and a — 2, a rate 

547 of ^450 triggered events/year is expected, which 

548 reduces to ^30 events/year for a = 1. The energy 



0.5 - 

n r i i . i • i ' • i . , 

17.5 18 18.5 19 19.5 20 

log u (Energy [eV]) 

Figure 16: Expected number of triggered events per year 
as a function of energy, from a realistic simulation of the 
MIDAS telescope. A microwave flux density F Ie { = F® , and 
a coherence parameter a = 1 were used to parameterize the 
EAS microwave emission in the simulation. 

549 spectrum of the triggered events is shown in Fig. [16] 

550 for the latter case. 

551 6. Data taking performance 

552 The MIDAS telescope underwent an extensive 

553 period of commissioning during several months in 

554 2011 at the University of Chicago, which provided a 

555 validation of the overall design and a test of the per- 

556 formance and duty cycle of the detector in a partic- 

557 ularly challenging environment for RF interference, 
ssa As a matter of fact, trigger rates were found to be 

559 significantly higher than those expected from ran- 

560 dom fluctuations. The SLT rate due to accidental 

561 triggers, rbkg, is estimated to be 0.3 mHz: 

562 r b k g = Np a tt • n pix • (r FLT ) npix (T) npi *~ , (12) 

563 where N pat t = 767 is the number of SLT patterns, 

564 n p ; x = 53 is the number of pixels in the MIDAS 

565 camera, tflt = 100 Hz is the pixel FLT rate, and 
see t = 10 fj,s is the coincidence time window. 

567 The background rate of SLT events during data 
sea taking was well above the estimate of Eq. [12] and 

569 highly variable, ranging from 0.01 Hz to 2 kHz. 

570 The major source of background was found to 

571 originate from airplanes passing over the antenna 

572 on their way to a close- by airport. Radar altimeters 

573 on board of these aircrafts operate just above the C- 

574 Band frequency, and, while suppressed by the MI- 

575 DAS receiver bandwidth, their emissions are strong 

576 enough to produce a sudden rise of the RF back- 

577 ground in many neighboring channels (or even in 
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Figure 17: Event display of a candidate detected by the MI- 
DAS telescope. See caption of Fig. |15| for a description of 
each panel. Although the timing characteristics of the event 
are compatible with those of a cosmic ray shower, correlated 
signals in off-track pixels do not allow for unambiguous iden- 
tification. 



the whole camera), with a corresponding increase 
of the first and second level trigger rates over sev- 

602 

eral tens of seconds. 

603 

Most of the remaining background was charac- 604 

terized by very fast (<200 ns) transient RF sig- 605 

nals, which trigger simultaneously from a few pix- 606 

els up to the entire camera. The frequency of these 607 

events varied greatly during the data taking period, 608 

with no particular correlation with day/night or 609 

week/ weekend cycles. While their exact origin was 610 

not clear, they were likely to have an anthropogenic 611 

origin in the vicinity of the detector. 612 

In the early stage of commissioning, a simple 6u 

veto, which inhibited the trigger system whenever m 

the SLT rate exceeded a given value, was imple- m 

mcnted in the FADC firmware to avoid flooding the eie 

DAQ with background events. Under this running 6i? 

conditions, the average live time of the detector was oia 

60%. A significant improvement was achieved with ei9 

the installation of a bandpass filter designed to re- 020 

ject the radar frequency band, which increased the 021 




Figure 18: Event display of a 4-pixel candidate detected by 
the MIDAS telescope. See caption of Fig. |15| for a description 
of each panel. Several candidates with similar characteristics 
- small signals and short tracks - were found in the data. 



effective duty cycle of the MIDAS telescope to bet- 
ter than 95%. 

The MIDAS telescope experienced very different 
weather conditions in Chicago - from high temper- 
atures during summer, to heavy rain and storms, 
to snow and ice. Detector operation and data tak- 
ing, as well as the telescope sensitivity, were found 
to be remarkably stable during the entire period of 
commissioning. 

The main objectives of the MIDAS telescope at 
the University of Chicago - to validate the telescope 
design, and to demonstrate of a large detector duty 
cycle - were successfully accomplished. In addition, 
a first search for events compatible with EAS was 
performed. Examples of candidates are shown in 
Figs. 



17 and 18 



However, the characteristics of 
these events, which have signals close to the trig- 
ger threshold, are also compatible with a tail of the 
overwhelming background noise. No strong con- 
clusion can be drawn on the origin of these events 
until a coincident detection with well-established 
techniques - for example at the Pierre Auger Ob- 
servatory - is performed. A detailed account of the 
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622 search and corresponding limits on microwave emis- 

623 sion from EAS are given in |28) . 

624 7. Summary and conclusions 

625 The MIDAS telescope - a prototype of an imag- 

626 ing detector for microwave emission from EAS - has 

627 been built and successfully operated at the Univer- 
62a sity of Chicago. The telescope's design is based 

629 on inexpensive off-the-shelf microwave components, 

630 and a custom-made digital electronics and trigger 

631 system. The absolute calibration of the detector 

632 was established with measurements of the Sun and 

633 other astrophysical objects as calibrated sources of 

634 microwave radiation. The sensitivity over the fo- 

635 cal plane was determined with an EM simulation 

636 of the telescope validated by measurements of the 

637 Sun transit in the telescope's field of view. The 

638 sensitivity of the MIDAS detector to EAS has been 

639 studied with Monte Carlo simulations which include 

640 a realistic parameterization of the detector based 

641 on these calibration measurements. Several tens to 

642 several hundreds of events per year are expected for 

643 microwave flux intensities as suggested by labora- 

644 tory measurements |llj . 

645 Several months of data taking in Chicago demon- 

646 strated that MIDAS can be reliably operated with 

647 minimal maintenance, and reach a duty cycle close 

648 to 100% even in an environment with high levels of 

649 RF background noise. Performances are expected 

650 to further improve in the radio quiet environment of 

651 the Pierre Auger Observatory, where MIDAS mea- 

652 surements in coincidence with the FD and SD will 

653 be essential to demonstrate the potential of this 

654 novel detection technique of UHECRs. 
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